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ABSTRACT 


Context. The Galactic center is the closest region in which we can study star formation under extreme physical conditions 
like those in high-redshift galaxies. 

Aims. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. 
Methods. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3 — 2 para-formaldehyde (P-H 2 CO) 
transitions. We used the 32, 1 — 22 , 0 / 3o,3 — 2o,2 line ratio to determine the gas temperature in n ~ 10^ — 10^ cm“ 
gas. We have produced temperature maps and cubes with 30^^ and 1 km s“^ resolution and published all data in FITS 
form*. 

Results. Dense gas temperatures in the Galactic center range from ~ 60 K to > 100 K in selected regions. The highest 
gas temperatures Tq > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km s~^ 
and 50 km s“^ clouds, and in “The Brick” (GO.253+0.016). We infer an upper limit on the cosmic ray ionization rate 
CcR < 10“^^ S“L 

Conclusions. The dense molecular gas temperature of the region around our Galactic center is similar to values found in 
the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the 
dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed 
temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, 
so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed 
to be high in the inner ~ 75 pc are confirmed to be high in the entire CMZ. 


1. Introduction 


The central region of our Galaxy, the central molecular zone 
(CMZ), is the nearest location in which we can study star 
formation in an extreme environment with pressure, tur¬ 
bulent Mach number, and gas t emperature much highe r 
than in local star-forming regions ( Morris fc Serabyn|1996 ). 
While there have been great leaps in our understanding of 
star formation in local molecular clouds in the past decade, 
there remain many unanswered questions about how star 
formation changes as gas becomes denser, more opaque, and 


* The data can be accessed from doi: 10.7910/DVN/27601 and 
are available from CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/ 
qcat?J/A+A/ 


more turbulent, as it mos t likely was in galaxies at the p eak 
of cosmic star formation (Kruijssen 


m galaxies at the p e 

fc Longmore||2Q13| ). 


The Galactic CMZ has attracted a great deal of atten¬ 
tion recently because it has a much higher ratio of dense 
gas ma ss to star formation tracers than the rest of the 
galax y ( Longmore et al.|2013bl 20T^ 2013al Kruijssen et al. 


2014 Yusef-Zadeh et al. 2QQ9[ ~ Immer et al. 2Q12| ). TTe 


star formation rate in this region therefore appears to be 
suppressed relative to the expectations from nearby star¬ 
forming regions and nearby galaxy disks, where the star for- 


sity ( 

Kennicutt|1998 Kennicutt & Evans|2012 

Leroy et al. 

2013 

Heiderman et al.|2010|). Within the CMZ 

, the central 
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Tlu 


^ 100 p 

SQfue[1995 
gas. Longmore 


ring” including the dust ridge (Lis et al. 1991 


Molinari et al.|20TT ) contains most of the dense 
et al.j ( |20i3bp recently proposed that the 


clouds along the “ring’ (which is m ore accurately described 


as a stream; 


In principle, this would allow us to study the time evolution 
of protoclust er clumps given an or bital model to describe 
their motion ( Krnijssen et al.|[2015 ). 

The thermal properties of this gas are a crucial compo¬ 
nent for understanding the conditions governing both star 
formation and interstellar chemistry. The gas temperature 
sets the sound speed in the gas and therefore the Mach 
number within a turbulent medium. It also governs the 
thermal Jeans mass. Regardless of whether turbulence or 
gravity controls the fragmentation scale, the gas tempera¬ 
ture is clearly important. 

Prior gas temperature measurements in the CMZ pri¬ 
marily used the popular ammonia metast able inversion 
transition thermometer (NH 3 (1,1) - (7, 7); Glisten et al. 
1981[ [Manersberger et al.|[l986| [Hiittemeister et al.||1993[ 
Ott et al. 2014), which is sens itive to mode rate density 


gas (^n(^H 2 


10 '=' - 10 ^ 


,-3. 


Shirley 


et al.| 

|2014|), so confirmation of the hi 

gh temperatures seen 
used to assess the sig- 

in NiJ 
nifica: 

[3 (Mills & Morris|2013j) can be 

nee of that mechanism. [Ao et al. 

(2013) used the para- 


& Wootten 1993), which is sensitive to denser (n 
— ^ -- 


cm and hotter {Tq > 20 K, usually Tq 


10 ^- 10 ^ 
60 K) gas, to 

map out the inner ^ 75 pc of the CMZ and found high tem¬ 
peratures comparable to those found in previous studies. 

There is a long-standing problem that the observed gas 
and dust temperatnres do not agree throughout the CMZ 


2011 


Ao et al. 2013 


hoton-dominated re- 


dCiisten et al.| 19811 Molinari et al. [2 

Clark et al ||2013[ |Ott et al.[|2014p . P 
gion TFonr models predict tnat gas and dust temperatures 
should ma tch at high column densities (Hollenbach & Tie- 
lens 1999). The observed discrepancy represents a signifi¬ 


cant problem for understanding the importance of gas tem¬ 
perature and pressure in the evolution of molecular clouds. 
The dust temperature is readily estimated via multiband 
continuum mapping and therefore is often used as a proxy 
for the dense gas temperature. While the discrepancy is not 
theoretically a problem, as the gas and dust may remain col- 
lisionally uncoupled up to moderately high densities given 
a high enough heating rate, it is an important empirical dif¬ 
ference between CMZ and Galactic disk clouds. Dust and 
NH 3 derived temperatures are usually a ssumed and often 
observed to agree in Galactic disk clouds ( Pillai et al.|2QQ6 


Dunham et al. 2010 Jnvela et al.||201^ Battersby et al. 


Ass uming the lAU-recommended d = 8.5 kpc (|Ghez et al. 


2008 


2013). 


IGillessen et al.||2009| |Reid et al.| 12009] [Gillessen et al. 


Krnijssen et al.|2015 ) represent a time-ordered 
path from starless to star-forming clouds, from “The Brick” 
(GO.253+0.016, cloud a) through the lettered dust ridge 
clouds (clouds b,c,d,e,f Lis et al. 1999), ending at Sgr B2. 


2014), but the difference can no longer be ignored in the 
Galactic center. 

Within both our own Galactic center and many nearby 
galactic nuclei, there is controvers y over which he ating 
mechanisms dominate. In the CMZ, Ao et al. (2013|) were 
able to rule out X-ray heating as an energetically important 
mechanism, leaving cosmic rays and turbulence (mechani¬ 
cal heating) as viable candidates. Studies of the CO spec¬ 
tral line energy distribution (SLED) around Sgr A* and 
Sgr B 2 arg ued for mechanical heating and UV heating, re - 


spectively (Goicoechea et al. [2013 [Etxaluze et al 


Nearby infrared-lumi nous galaxiesnave properties th at ro 


2013). 


semble these regions (Kamenetzky et al. 2012, 2014), but 


the debate about the dominant heating mechanism remains 
open in many galaxies, allowing for photon-dominated, X- 
ray- dominated, mechanical/tnrbnlent, or cosmic ray heat¬ 
ing (Loenen et a l. |2008 van der Werf et al. 2010[ |Pa- 


padopoulos et al.||2Qll[ |Meijerink et al.||2011[ |Bayet et al. 

2011[ Mangum et al.||2013 Papadopoulos & Thi 2013).^ 

We have expanded the H 2 CO mapping project of |Ao| 

et al. (2013) to a ^ 5 x larger area, including the dust ridge 


2015). While this 


thermometer is generally reliable in cold regions (e.g., the 
(1-1)/(2-2) ratio is accurate to Tq < 40 K), the population 
of the higher energy states may be affected by formation 
heating, a mechanism recently discovered to affect H 3 O+ in 
the CMZ and suggested to affect NH 3 , in which the excess 
energy from molecule formation leaves a significant frac- 


and the positive-longitude turbulent clouds thought to be 
associate d with the intersection point between the xl and 
x2 or bits ( [Rodriguez-Eernandez et al.|2006 Riqnelme ^ al. 
' 2013[ ). We describe the new observations in Section In 
Section we describe the analysis process used to extract 
signal from the data cubes. In Sectionj^ we describe how we 
derive temperatures from the line ratios. We discuss various 
implications of our data in Section In the appendices, 
we describe further details of the data reduction process, 
include additional tables, and provide the source code for all 
aspects of this project, from data reduction through figure 
generation for the paper. 


2. Observations data reduction 

2.1. Observations 

We observed the GMZ from —0.4° < ^ < 1.6° with the 


APEX -1 (SHEI) instrument (Vassilev et al. 2008) on the At- 
acama Path finder Experiment (APEX) telescope (Glisten 
et al. 2006[ ) using the extended bandwidth East Eonrier 
Trans form Spectrometer (XEETS) backend (Klein et al. 
2012). The observations were performed in service mode 
and were spread out over two years. The time was divided 
into 25 hours in June 2013, 75 hours in April-July 2014, 
and 50 hours in October 2014. A final set of observations 
was taken in April 2015; these data are not included in 
the analysis in this paper but are provided in the delivered 
EITS files. The time was split between the ESO (E-093.G- 
0144A; 50h, E-095.G-0242A, 25h), MPIfR (M-091.E-0019 
and M-093.E-0009; 75h), and OSO (E-093.C-0144A; 25h) 
queues. 

The 2013 observations were taken in 4' x 4' patches, 
and the frequency range covered was 217.5-220 GHz and 
216-218.5 GHz in the two spectral windows. Scans were 
performed along lines of constant RA and Dec. These ob¬ 
servations used the same observing strategy and off position 


Ao et al 


The 2014 observing strategy was modified to use larger 
8' X 8' scans along lines of Galactic latitude and longitude. 
The frequency range was also shifted to cover windows over 
217-219.5 and 218.4-220.9 GHz, thus including the bright 
and G^^O 2-1 lines. Three off positions were used 


Article number, page 2 of 41 
















































































































































































Ginsburg et al: APEX CMZ H2CO 


for these observations: (a) 17:52:06.854 -28:30:31.32, (b) 
17:43:53.890 -28:07:04.68, and (c) 17:48:11.934-29:44:41.83; 
position (b) exhibits some emission from local clouds 

near 0 km s“^ but no other emission lines, while both (a) 
and (c) appear to be entirely clean. 

Additionally, for the p-H^C O data, we incorp orated the 
41 hours of observations from Ao et al.| (|2013|) using the 
older fast fourier transform spectrometer (EFTS) backend. 
These data covered 2 GHz of bandwidth, including all three 
of the P-H 2 CO lines and SiO 5-4, but they did not cover 
the CO isotopologue lines. 

The detected lines were P-H 2 CO 3 o ,3 — 2 o ,2 218.22219 
GHz, P-H 2 CO 32 , 2 - 22,1 218.47563 GHz,’p-H 2 CO 32 ,i- 22 ,o 
218.76007 GHz, SiO 5-4 217.10498 GHz, CH 3 OH 42,2 - 3 i ’2 
218.44005 GHz, OCS 18-17 218.90336 GHz, HC 3 N 24-23 
218.32471 GHz, SO 65-64 219.94944 GHz, HNCO 10o,io- 
9o,9 219.79828 GHz, 2-1 219.56036 GHz, and i^c’o 2- 
1 220.39868 GHz. All listed frequencies are rest frequencies, 
but we do not include many, e.g., CH 3 CN 12-11, that are 
in the observed band but were detected only in Sgr B 2 , 
since complete surveys of the rich m olecular heimat in this 
band have been com piled elsewhere (Nummelin et al. |1998i 
Belloche et al.|2013 ). Of the detected lines, the SO, HNCIt 


C^^O, and ^^CO were not covered by the 
data but were in the newer XFFTS data. 


Ao et al. 
WEIF we 


2013) 


o not 


analyze these lines in this paper, we provide data cubes that 
include them. 

The raw data were acquired with 32768 spectral chan¬ 
nels in each window, yielding 0.1 km s“^ resolution. We 
expect to see no lines narrower than a few km s“^ in the 
CMZ, particularly not with the relatively shallow observa¬ 
tions we have acquired. We therefore downsampled the data 
by a factor of 8 to 0.8 km s“^ resolution prior to resampling 
onto a 1 km s“^ grid to make the data more manageable. 

< 


The system temperature ranged from 120 < T, 


sys 


200 K for the majority of the observations, with a mean 
Tsys = 165 K. A small fraction (10%) of observations were 
in the range 200 < Tsys < 300 K. There were also a very 
small number (< 1 %) with much higher system temper¬ 
atures (300 < Tsys < 750 K) that were flagged out (see 
Section 2 . 2 . 2 | ). 


2.2. Reduction 
2.2.1. Calibration 

Calibration was performed at the telescope using the stan¬ 
dard APEX calibration tools These yield flux-calibrated 
spectra at each position with appropriate pointing informa¬ 
tion. Typical flux calibration uncertainties are ^ 10% and 
pointing errors a < 2 ". 

There was a significant calibration error discovered at 
the APEX telescope during a large segment of the 2014 
observing campaign (http://www.apex-telescope.org/ 
heterodyne/shfi/calibrat ion/calf actor/). The cali¬ 
bration errors were of order ~ 15 — 25%. We corrected 
the data using the correction factors suggested on the cal¬ 
ibration page, but the remaining calibration uncertainty 


^ See the observing manuals: www.apex-telescope. 

org/docmnents/public/APEX-MPI-MAN-0012.png 
WWW. apex-telescope.org/documents/public/ 

APEX-MPI-MAN-0013.png 


is higher in this data, ^ 15% total rather than the usual 

- 10%. 


2.2.2. Flagging bad spectra 

Spectra were removed if they showed excessive noise com¬ 
pared to the theoretical e xpectation given the measured 


system temperature. As in |Ao et ah (2013), the threshold 


was set to 1.5 x the theoretical noise from averaging two 
polarizations, i.e. STsysi^^^J^texp) 


-0.5 


where texp is the ex¬ 


posure time (integration time) per spectrum in seconds and 
Au is the channel width in Hz. This approach resulted in 
^ 0 . 2 % of the data being removed. 

At the position of Sgr B 2 , the noise was significantly 
higher due to signal from the continuum source. We there¬ 
fore disabled this flagging in a 2.5' box around Sgr B 2 . 


2.2.3. Mapmaking 

The maps were made by computing an output grid in Galac¬ 
tic coordinates with 7.2" pixels and adding each spectrum 
to the cube at the appropriate pixe|^ location. In order to 
avoid empty pixels and maximize the signal-to-noise, the 
spectra were added to the grid with a weight set from a 
Gaussian with FWHM = 10", which effectively smooths 
the output images from FWHM ^ 28" to 30". See 
Mangum et al. (2007) for more detail on the on-the-fly 
mapping technique used here. The spectra were averaged 
with inverse-variance weighting. The full weight equation 
is therefore: 


Wi = exp 


- {{x - Xif + {y- ViY 


( 1 ) 


where Xi and yi are the coordinate of the pointing in pixel 
space, apix ~ 10"/v^8To^/(7.2"pix“^) = 1.38, and arms is 
the standard deviation measured along the spectrum in a 
signal-free frequency region. 

The position-position-velocity (PPV) cubes were cre¬ 
ated with units of brightness temperature on the corrected 
antenna temperature scale (T^h The main beam efficiency 
is rjmb = 0.75 (gain ^ 39 Jy/K)r These values are noted in 
the FITS headers of the releas^ data. 

The maps achieved a depth of a = 50 — 80 mK in 1 
km s“^ channels and FWHM = 30" beams. The noise is 
slightly lower (about 15%) than expected from the online 
APEX calculator because our noise measurements are made 
in moderately smoothed maps. 


2.3. Baselining 

The data showed significant baseline structure, leading to 
large-scale correlated components in the resulting spectra. 
The baselines were removed over a velocity range —150 < 
'^LSR < 150 km s~^ by the following procedure: 

1 . identify bright PPV regions in the P-H 2 CO 3 o ,3 — 2 o ,2 
line using the technique described in Section 

^ We use the term ‘pixel’ to refer to a square data element 
projected on the sky with axes in Galactic coordinates. The 
term ‘voxel’ is used to indicate a cubic data element, with two 
axes in galactic coordinates and a third in frequency or velocity 
^ http://www.apex-telescope.org/telescope/efficiency/ 
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Fig. 1. Spectra toward one H 2 CO-bright cloud, the 20 km s“^ cloud {£ = 359.892, b = —0.074), over the full APEX band. All 
detected lines are identified. The feature overlapping the 2-1 line at about 219.5 GHz is the line from the upper 

^X)f|4?P There are artifacts at the band edges, e.g., below 217.0 GHz, because these spectra average a 
few different tunings withmre accounting for missing data at the band edges. 






































Ginsburg et al: APEX CMZ H2CO 


2. mask the regions identified as 80,3 — 2o,2-bright for each 
other spectral line data cube 

3 . fit a polynomial to the un-masked data 

4 . subtract the fitted polynomial 

We used a 7 th-order polynomial for this process because 
lower order baseline removal left a significant and patterned 
residual. A more detailed examination of the baseline re¬ 
moval process is described in Appendix [A| 


becomes undefined. We therefore limited our ratio measure¬ 
ments to regions with significance > 5cr in the 80,3 — 2o,2 
line and used that line as the denominator in the ratio. By 
using the brighter line as the denominator, we avoid divide- 
by-zero numerical errors. 


3.2. Region extraction 


3. Signal extraction masking 


We use the method descr ibed partially in Ao et al. ( 2013 ) 
and more thoroughly in Dame ( 2011 ) to mask the data 
cubes at locations of signihcant signal in the brightest line. 
A noise map was created by computing the sample stan¬ 
dard deviation over a 200 km s“^ range in which no signal 
was present. We use the P-H2CO 80,3 — 2o,2 line, which is 
the brightest of the H2CO lines, to create the mask by the 
following procedure: 


1 . Smooth the data with a Gaussian of width two pixels in 
each direction (spatial, a = 14 . 4 ", and spectral, a = 2 
km s“^) 

2 . Create an inclusive mask of all pixels with brightness 
Ta > 2(7 in the smoothed map 

3 . Grow the mask from the previous step by one pixel in 
each direction (this is known in image processing as bi¬ 
nary dilation) 


The P-H2CO 3 o,3 — 2 o,2 mask was then applied to the 

32.1 ~ 22,0 aiid 82,2 — 22,1 cubes created with the same PPV 
gridding. There is some overlap between the CH3OH 42,2 — 

31.2 line and the P-H2CO 82,2 — 22,1 line in PPV space, so 
we shifted the P-H2CO 80,3 — 2o,2 mask to the velocity of the 
CH3OH line in the P-H2CO 82,2 — 22,1 cube and excluded 
those pixels expected to contain CH3OH emission. However, 
inspection of the 82,2 — 22,1 cube revealed that significant 
chunks of the 82,2 — 22,1 signal were excised by this masking, 
so all further analysis excludes the 82,2 — 22,1 line. 

Integrated intensity maps are shown in Figure for 
the 3 o,3 — 2o, 2 map masked with the above method and in 
Figure for the 82,1 — 22,0 map using the mask created 
from the 80,3 — 2o,2 data. 


3.1. Ratio maps 

We have created P-H2CO S'zy(32,i —22,0)/S'zy (80,3—20,2) ratio 
maps as a first step toward building temperature maps. In 
general, a higher P-H2CO 3 ^( 32 ,1 - 22,o)/5'zy(3o,3 - 20,2) in¬ 
dicates a higher temperature, and therefore the ratio maps 
are a good proxy for relative temperature. We displa y four 
different varieties of ratio map in Figures [ 3 | and |E.l] 

Figure shows maps produced by creating ratio cubes, 
in which each voxel has the value Ri = P-H2CO 5'zy(32,i — 
22,o)/5'zy(3o,3 — 20,2), then averaging across velocity with 
constant weight for each voxel in Figure and with the 
weight set by the 80,3 — 2o,2 brightness value in Figure!^. 
In both cases, the mask described in Section was applied 
to the cubes before averaging over velocity. 

Ratio maps are subject to dramatic uncertainties in the 
low signal-to-noise regime, which is inevitably reached in 
some parts of the map. The error distribution on the ratio 
becomes non-Gaussian as signal approaches zero, approach¬ 
ing a pathological Cauchy distribution in which the mean 


In order to extract higher signal-to-noise measurements on 
selected regions, we broke the data down into subsets using 
both a by-eye region selection drawing regions using ds 9 


and a more systematic ap proach using a dendrogr am-based 
clump finding algorithm ( [Rosolowsky et ar]| 20 Q 8 , http:// 
dendrograms. org/). The dendrogram clump hnding uses a 
watershed decomposition like other clump finding methods, 
but it also extracts clumps at different brightness scales 
(different heights in the watershed) and tracks how nested 
clumps are related to each other. Within a hierarchy of 
nested clumps, the brightest and smallest is a deaf’ and all 
associated larger structures are ‘ancestors’ (roots or trunks 
and branches). We refer to the clumps identified in this 
manner as ‘dendrog ram-extracted’ clumps fr om here on. 


Figures 1 and 8 in Rosolowsky et al. ( 2008 ) give visual 


representations of this hierarchy. 


We ran the dendrogram extraction analysis on the 

3o,3 — 2o,2 PPV cube data. We assumed a uniform noise 
per pixel (Tt = (ctt) = 0-07 K, which corresponds to the 
mean noise per field. We enforced thresholds of a mini¬ 
mum number of voxels = 100, a peak threshold of 

/ = 8(7, and a splitting threshold A = 2cr, which defines 
how large the brightness difference between the peak pixel 
and an adjacent pixel must be in order to divide the clump 
into two sources. The exact values of these parameters are 
not particularly important, as we are interested in general 
trends with size-scale and galactocentric distance, but we 
caution against overinterpretation of the resulting catalog 
as the number of sources and their size and distribution 
can change dramatically with small changes in the selected 
parameters. The number and location of voxels included 
in the catalog is, however, relatively robust against small 
parameter changes. 


For each extracted clump, we measured the correspond¬ 
ing integrated and peak emission in the 80,3—2o,2, 82,1—22,0, 
and cubes. We also extracted the mean dust tem- 


(Molinari et al. 

2010 Traficante 

s were performec 

on background- 


snbtracted data using a n approach originally described in 
Battersby et al.| (2011|), which includes details of the fit- 
ting process and specification of assumed physical param¬ 
eters. We also extracted dust column density and tem¬ 
perature using a more naive pixel-by-pixel approach with 
no background subtraction (http://hi-gal-sed-f itter. 
readthedocs.org/en/latest/) and found the two to be 
consistent. We assumed a dust to gas ratio o f 100 and 
dust opacity = 4 g cm“^, the same as Battersby 


et al. ( 2011 ) and equivalent to that assumed in the Bolo- 


P — 2 (^271 

AC.Hz = 1-1 g cm 2 

Bally et al. 

2010 Aguirre 

et al.pOll 

Ginsburg et al .||2013 

). 
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Fig. 2. In teg rated intensity (moment 0) maps of the (a) P-H 2 CO 80,3 — 2 o ,2 and (b) 82,1 — 22,0 lines. The cubes were masked prior to integration using the method described 
in SectionSelected region names are labeled. The thin black contours are at a column iV(H 2 ) = 5 x 10^^ cm“^ from the Herschel SED fit maps. 
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cubes. The Sgr B2 peaks exhibit H 2 CO absorption and therefore do not contain reliable measurements. Higher ratios (red) correspond to higher temperatures, 
(b) The same as (a), but the average along velocity has been weighted by the 80,3 — 2 o ,2 brightness 
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3.2.1. Manual region extraction 

The by-eye selected regions were extracted less systemati¬ 
cally, instead focusing on peaks in the 80,3 — 2 o ,2 or 82 , 1 — 22,0 
emission. These selections were made in two dimension 
rather than in the full PPV space. We included a series 
of large square 8 ' x 8 ' apertures corresponding to each ob¬ 
served field. For each region extracted in this fashion, we fit 
a 1-8 component model as described in Section [4^ where 
each component consists of a 6 -parameter multi-Gaussian 
model. An example spectrum extracted from a small re¬ 
gion around “cloud c” is shown in Figure The extracted 
regions are shown in Figure [BT| 


4. Analysis: temperature determination 

In this Section, we describe the techniques used to derive 
gas temperatures. In Section |4.1| we model the line ratios 
with a radiative transfer code. We fit Gaussians to spectra 
extracted from the cubes and use these fits to measure line 
ratios in Section [T^ In Section [T^ we apply the same mod¬ 
eling a ppro ach to the dendrogram-extracted clumps f rom 
Sect ion | 8 . 2 [ We use the modeling results from Section |4.2| 
and |4.8| tojustify a simplified modeling approach that we 
apply to the entire map of the CMZ in Section |4.4[ 


4.1. Radiative transfer modeling 


We use RADEX (van der Tak et al.||2007|) and a related 
solver (Fujun Du’s myRadex; https : //github. com/f j du/ 
myRadex) to create model grids for the P-H 2 CO lines over a 


grid of 20 densities (n = 10 ^‘^- 10 ^ cm ^), a fixed assumec 


line gradient of 5 km s“^pc“^, 80 H 2 CO column densities 
A^(H 2 C 0 ) = 1 X 10^^ — 1 X 10^^’^ cm“^, and 50 temper¬ 
atures from 10 to 850 K. The grid was created using the 
pyradex wrap per of RADE}([^ and collisi on rates retrieved 
from LAMDA ( Wiesenfeld fc Fanre|2018 ). We then upsam¬ 
pled these grids by spline interpolating between grid points 
to acquire a high-resolution {Tq xnx N = 250 x 100 x 101) 
grid covering 10 < < 850 K, 10^-^ < n < 10^ cm“^, and 

10^^ < N < 10^^’^ cm“^ km s“^pc“^. Eor our analysis, we 
used the myRadex grid; it is consistent with the RADEX grid 
over most of parameter space but suffers fewer numerical 
instabilities over the region we are interested in. 


4.2. Spectral modeling 


We extracted spectra averaged over selected regions and fit¬ 
ted a 6 -parameter model to the full 218-219 GHz sp ectral 
range using pyspeckit (jGinsburg & Mirocha| 2011). The 
fitted parameters are the amplitude of the H2GO 80,3 — 2o,2 
line, the centroid velocit}!^ the line width (cr), the ratio 
Ri = P-H 2 CO 5 'z.( 82 ,i - 22 ,o)/S'iy( 8 o ,3 - 20 , 2 ), the ratio 
R2 = P-H 2 CO 5 'zy(82,2-22,i)/5'iy(82,1-22,0), and the ampli- 


^ The line gradient is selected because observed line widths have 
FWHM larger than 5 km s“^ in 30'' (1.2 pc) beams throughout 
the CMZ. The temperature measurements are robust against 
this assumption, with a typical change of dTjT < 25% for a 
gradient 4x larger or smaller. 

® https://github.com/adamginsburg/pyradex 
^ We assume the CH 3 OH and H 2 CO lines have the same cen¬ 
troid velocity, and we find no counterexamples to that assump¬ 
tion within our data. 
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tude of the CH 3 OH 42,2 —81,2 lin^ where Sj^ is the surface 
brightness. For these spectra, we performed a different base¬ 
line subtraction to that described in Section 12.81 We first 
fit the described 6 -parameter model, then fit a spline curve 
to the residual smoothed to a 50-channel (50 km s“^) scale. 
We then re-fit the model to the data. In all but the lowest 
signal-to-noise cases, this approach resulted in consistent 
ratio measurements between the two fits and with the other 
methods. Figure shows two example spectra that exhibit 
significantly different ratios with the spline baseline fit in¬ 
dicated. We use the baseline-subtracted measurements for 
further analysis. 


For each extracted region (see Section 8.2), we created 


grids over the modeled temperature-density-column den¬ 
sity parameter space. We measured with five indepen¬ 
dent constraints: the line ratio Ri, the H 2 CO abundance, 
the total column density of H 2 , a lower limit on the H 2 
volume density, and an upper limit on the beam filling 
factoij^ We use the line ratio as the primary constraint 
rather than line brightnesses to avoid uncertainties due to 
the filling factor of the emitting gas: even for diffuse clouds, 
the filling factor of the emitting regions may be // << 1 
if the emission is isolated to compact shocked regions, as 
expected if highly supersonic AI > 10 turbulence is en¬ 
ergetically dominant in the clouds. We use an abundance 
A^(p-H 2 CO)/A^(H 2 ) = X(p-H 2 CO) = allowing 


H 2 CO abundance dGinsburg et al.| 

12018 

Carey et al.||1998t 

Wootten et al. 

1978| Mundy et al.j 

00 

). To constrain the 

lel dust maps (see 

total column c 

ensity, we use the 1 

lerscl 


Section 8 . 2 ) to derive an H 2 column density, which has a 
nominal ^ 2 — 8 x uncertaint}!^ We treat the error as lOx 
(instead of ^ 2 — 8 x) to account for the additional un¬ 
certainty in the abundance of P-H 2 CO. The effective range 
of A^(p-H 2 CO) is therefore < N < cm“^, al¬ 
lowing for < A^(H 2 ) < cm“^ and Tier in log 

abundance. To constrain the dust density, we assume the 
selected area has a mass given by the Herschel dust-derived 
mass and a spherical volume, which sets a conservative (but 
nonetheless highly constraining) lower limit on the volume 
density. 

These constraints are shown projected onto the three 
planes of our fitted parameter space in the multi-paneled 
Figures rD.l||D.4| for one example region. The fitted param¬ 
eters are displayed in Figure [P^ 


While R 2 can provide significant constraints on various 
parameters, we do not use it for our general analysis because 
of the ambiguity imposed by the overlap with the CH3OH 
42,2 — 81,2 line. Blending between the CH3OH 42,2 — 81,2 
and H2CO 82,2 — 22,1 line occurs frequently, in ~ 1/4 of all 
sightlines. 


^ The parameters R 2 and the CH 3 OH brightness are essential 
for fully fitting the spectrum, but they are not used for further 
analysis. 

^ Where the model predicts a brightness temperature higher 
than observed, we can invoke a low filling factor to explain the 
discrepancy. If the model predicts a brightness temperature too 
low, we can rule it out. 

The statistical uncertainty is much smaller than the system¬ 
atic unknowns, including the dust opacity per unit mass, the 
presence of multiple line-of-sight cloud components, and the 
presence of multiple temperature components within any given 
cloud. 
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Brick SW 



ds:G0.42+0.03box 



(b) 


Fig. 4. Simultaneous fits to four lines, from low to high frequency, P-H 2 CO 80,3 — 2 o, 2 , CH 3 OH 42,2 — 81 , 2 , P-H 2 CO 82,2 — 22 , 1 , and 
P-H 2 CO 82,1 - 22,0 for (a) The Brick {i = 0.241, b = 0.006, r = 25")and (b) a box centered on 00.42+0.08 {i = 0.415, b =’o.082, 
80'^ X 80'^). The fitted parameters shown in the legend are the amplitude of the H 2 CO 80,3 — 2o ,2 line, the centroid velocity, the line 
width (cr, not FWHM), the ratio Ri = P-H 2 CO S'zy( 82 ,i — 22 ,o)/*S'iy( 8 o ,3 — 20 , 2 ), the ratio R 2 = P-H 2 CO 5 'zy( 82,2 — 22 , 1 )/iSiy( 82,1 — 22 , 0 ), 
and the amplitude of the CH 3 OH 42,2 — 81,2 line. These spectra have significantly different ratios and therefore derived temperatures: 
The Brick has Tg ^ 100 K while GO.42+0.08 has Tg ~ 40 — 50 K. The bottom black spectrum is the residual, and the orange 
wiggly spectrum shows the spline fit used to remove baseline ripples. The OGS 18-17 218.90886 GHz and HG 3 N 24-28 218.82471 
GHz lines are also detected, but not fitted, in The Brick. In these examples, we are showing cases where the GH 3 OH 42,2 — 81,2 
and P-H 2 GO 82,2 — 22,1 lines are not too badly blended, but along many lines of sight, especially those within ~ 20' of Sgr B2, 
they are irretrievably blended. 


c:G0.38+0.04 



Fig. 5. Fitted spectrum of “cloud c". 


The fitted parameters and their corresponding errors are shown in the legend as in Figure]^ 


4.3. Clump modeling 


We derived temperatures for each dendrogram catalog 
clum p using the same measured quantities as in Sec tion 
4.2 as constraints. The described in Section 4.1 was 


used to extract the temperature and associated error. 


We plot the derived temperature for all clumps against 
the measured line ratio in Figure This figure includes 
both ‘leaves’ (the most compact of the nested clumps) and 
‘ancestors’ (any clump that contains a smaller clump within 
it) from the dendrogram extraction. The figure illustrates 
the effects of assuming higher or lower density on the tem¬ 
perature; the derived temperature is only weakly dependent 
on the assumed density. The lower density models are ex¬ 
cluded by the observed volume density lower limit (assum¬ 
ing spherical geometry) from the dust. We use an assumed 


density n(H 2 ) = 10^ cm ^ when a uniform assumed density 
is required in further analysis. 

4.4. Temperature maps 

Using the relations shown in Figure we converted the 
ratio map shown in Figurej^to a temperature map in Figure 
[^with a constant abundance X(H 2 CO) = 1.2 x 10“^ and 
two constant densities n = 10^ (Figure |^/b) and n = 10^ 
cm“^ (Figure]^). Because this map incudes the full line- 
of-sight integrated emission, there are many regions where 
multiple independent components are being mixed. 

These temperature maps do not have uniform reliabil¬ 
ity. There is some spatial variation in the noise, but the 
variation in the signal is much more important. In regions 
where the peak signal is low, the reliability of these maps is 
substantially reduced. We have also applied a mask to the 
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Ratio *S'(32 1—22 o)/‘S'( 3 o 3— 2 o 2) 

Fig. 6. The derived temperature vs the measured ratio P-H 2 CO Su{^ 2 ,i —‘^ 2 ,o)/Su{^o ,3 — ‘^ 0 , 2 ) for the dendrogram clumps. The big 
symbols represent leaves (compact clumps) and the small symbols represent larger and therefore lower-density ancestor structures. 
The temperature includes constraints from the assumed constant H 2 CO abundance and the varying column density. The points are 
color coded by signal-to-noise in the ratio with gray S/N < 5, blue 5 < S/N < 25, green 25 < S/N < 50, and red S/N > 50. 
The black lines show the modeled temperature as a function of Ri for three different assumed densities with a constant assumed 
abundance XH 2 C 0 = 1.2 x 10“^. While the measured values appear consistent with the n = 10^ cm“^ curve, the systematic errors 
permit any density in the range 10^ cm“^ < n < 10^ cm“^; this plot most importantly shows that n ~ 10^ cm“^ is not plausible 
and that the temperature does not differ very much for a given value of Ri over the plausible range of densities. 


lower panels of Figure to make the low-reliability regions 
blend into the background by adding a foreground gray 
layer with opacity set by the inverse S/N (see Appendix 

H- 

The effective range of detected temperatures is ~ 40 — 
150 K. While lower and higher temperatures appear in 
the map, those measurements should not be trusted. Be¬ 
low Tq ^ 40 K, all temperature measurements in the map 
are effectively Icr upper limit^^ This limit is due to a 
combination of sensitivity and excitation: below 40 K, the 
P-H 2 CO lines are observed to be faint; the 82,1 — 22,0 line 
is not detected. In regions of high column density, we de¬ 
tect both lines, and in all of these regions, we have mea¬ 
sured high temperatures {Tq > 40 K). Since Tq ^ 40 
K corresponds to Ri ^ 0.15, a 5cr temperature measure¬ 
ment would require Te( 32 ,i — 22 , 0 ) = (5 x 0.07K) = 0.35 
K and Tb( 3 o ,3 - 20 , 2 ) = (5 x 0.d7K)/0.15 = 2.33 K. This 
Tb{^o ,3 — 20 , 2 ) is slightly brighter than the brightest detec¬ 
tion in our map, and all regions detected with such high 
brightness temperatures also turned out to have higher 
Tb{^ 2 ,i — 22 , 0)5 nnd therefore higher Ri and temperature. 
The higher measured temperatures become lower limits 
above ^ 150 K because the P-H 2 CO line ratio is intrin¬ 
sically insen sitive to higher temperatures, as can be seen in 
Figure [P^ (rightmost column), in which the modeled line 
brightnesses are both nearly constant above Tq ^ 150 K. 


When averaged into regions in later sections, we measure some 
temperatures below this limit, albeit still at low significance. 
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5. Discussion 

5.1. Gas heating and thermal balance 


Following Ao et al. (2013) and Papadopoulos (2010), we 
examine the energy balance in the CMZ gas to understand 
the temperatures we measure. The two key processes they 
identified are cosmic ray and turbulent dissipation (shock) 
heating. We ignore all UV-driven heating processes in this 
section as we are interested in gas at high extinction {Ay > 
lOp^ Other processes, e.g., direct (proto)stellar heating 
or supernova bubble interaction may be relevant l ocally, 
but cannot heat the gas globally. |Ao et al^ ( |2013[ ) ruled 
out diffuse X-rays as a significant heating source, with an 
energy density ^ 3 orders of magnitude too low to explain 
observed temperatures. 

We determine the expected gas kinetic temperature by 
assuming that the heating and cooling processes are in 
equilibrium. This is generally a good assumption on the 
scales we observe, as typical gas cooling rates via line 
or dust cooling around Tq r\j 100 K and n ^ lO^cm ^ 
are dT/dt ^ 10~^ ^erg s“^cm“^ ^ 23(n/10^cm“^)K kyr“^ 
( Goldsmith||2001 and see below). The cooling timescale is 
shorter than the bimescales for any energy injection mecha¬ 
nisms that may plausibly heat the gas. The supernova rate 
is < 1 kyr“^ over the whole CMZ, and supernova rem¬ 
nants stop producing cosmic rays rapidly, so the super¬ 
novae are expected to contribute to a diffuse cosmic ray 
population rather than local stochastic heating incidents. 
Turbulence in these clouds is dyn amically driven on scales 
of 10-30 pc dKrnijssen et al.|2015 ). Assuming the turbulent 
dissipation scale is close to the sonic scale with a sound 
speed Cs = 0.76(Tg'/ 100K)^/^ km s“^, the gas temperature 
will be well mixed on milliparsec spatial scales on kiloyear 
timescales. 


niewicz 


Wing - fai Th i from work by the ProD i Mo team from |Le Bourlot 
et al. Wrathmall et al. (2007); Flower et al. (2000 )7 |Wol- 

zz et al. (|l998|. Input abundances were X{UO) = 10 “"^, 
X{C^^O) = 5 X 10 -^ X{C) = 5 X 10 -^ 
X{0) = 5 X 10“®, X{HD) = 10“^ and 


gradient (LVG) approximation for optical depth calcula¬ 
tion. The resulting curves for Tg{<Jv) are shown in Figure 
(see Appendix for further details). We have assumed 
a size scale L = 5 pc, a gradient dvjdr = 5 km s~^ pc“^, 
and a cosmic ray ionization rate (CRIR) (cr = 10“^^ s“^, 
with exceptions noted in the legend. We assume a dust 
temperature Td = 25 K, but we assume lower dust radia¬ 
tion temperatures because the dust is optically thin at long 
wavelengths in most regions. We have assumed a constant 
interstellar radiation field (ISRF) of Go, which should only 
weakly affect dust temperatures and has no effect on the 
gas at the high column densities being modeled. A higher 
ISRF may be appropriate but does not change our results, 
since a 1000 Go field provides the same gas heating rate 
as a (cR = 2.5 x 10“^^ s“^, a rate much lower than the 


measured CRIRs in the CMZ. IClark et al. (2013) similarly 
find that the ISRF has little effect on the gas temperature. 

A first key point is that we measure gas temperatures 
Tq ^ 60 K in many clouds everywhere throughout the 
CMZ. Since cosmic rays have very high penetration depths, 
we expect that cosmic ray heating is uniform everywhere 
within the CMZ, independent of column densit^f^ If this 
uniformity holds, then we can set an upper limit on the cos¬ 
mic ray ionization rate, (^cr ^ 10“^^ s“^, since a higher rate 
would result in a minimum gas temperature above what we 
have measured. The limit is consistent with many previ¬ 
ous measurements of the C RIR, though it rules out some 


of the more extreme values ( Yusef-Zadeh et ^|2013 Goto 


et al.|2013|). For a fiducial upper limit measurernent, we use 

- J A - Q,/, 1 - 1 ^ 


We expand upon the Ao et al. (2013) analysis by solving 
for the temperature assuming that both cosmic ray ioniza¬ 
tion and turbulent heating are present and that the dust 
temperature is nonzero. We assume the heating rate for 
cosmic rays is constant per cubic centimeter Tcr = nCcR^ 
and the heating rate from turbulence Tt^rb = ncr^IL, where 
a is the 3D velo city dispersion an d L is the length scale. We 
used DESPOTIC ( KrnmhoIz]|2Q14 ) to determine the equilib¬ 
rium temperature. We included CO, O, C+, C, 0 -H 2 , P-H 2 , 
and HDpq as line coolantj^ We used the large velocity 


While some emission likely comes from PDRs, we expect these 
regions to have small filling factors and typically lower densities 
than cloud centers. At densities n < 10^ cm“^, P-H 2 CO 3-2 
emission is very weak and would contribute little to the ob¬ 
served line emission. Volumes within the actively star-forming 
regions Sgr B2, Sgr C, and near the circumnuclear disk may be 
locally affected by strong UV heating from young and forming 
stars, but the bulk of the dense g as is not. Finally, PDR models 
(e.g., Hollenbach X Tielens|1999 ) predict that Td ^ Tp at high 
densities, while we show that Tg > Td', see Sectior5.3 

The collision rate files for H 2 and HD were provided by 


n = 10"* cm“^, dvjdr = 5 km pc“^, T^ad^D = 10 K, and 
CcR = 10“^^ s“^, which yields T^' = 72 K in the absence of 
any other heating mechanisms. This upper limit could be 
violated if the gas in all observed clumps is about an order 
of magnitude denser than the dust-derived average density, 
which is not likely. This limit leads to the conclusion that 
cosmic ray heating is either not dominant in the CMZ or 
not uniform, as uniform cosmic ray heating cannot satisfy 
the requirement for both high (> 100 K) and low (< 40 K) 
temperatures to be present in the CMZ. 

We have information on the linewidth, size scale, and 
density for each individual region, albeit with great uncer¬ 
tainties on each. To assess the role of heating by turbulent 
energy dissipation, we use these measurements to predict 
the gas temperature. We use the effective radius, which 
is the geometric mean between the moment-derived ma¬ 
jor and minor axee^ of the selected sources as the length 
scale, the measured linewidtlp^a^, the measured dust tem- 


Lis & Menten 
emission from 


V('"GO) = 5 x ir 

X(G+) = 5 X 10“^°, 
the ortho/para ratio of H 2 OPR = 3. 

We assume a very low abundance for C+, Xc+ = 5 x 10“^°, 
assuming that cloud shielding will prevent photoionization and 
collisions with H 2 wi ll destroy CR-gene r ated C + rapidly, which 
is confirmed by the Nelson X hanger (1999) chemical model. 


(11998 ) and Lis et al. (1999) observed bright C+ 
The Brick; we interpret this emission as coming 
from a heated outer skin of the cloud at lower densities not 
coincident with the observed H 2 CO emission. 

However, low-energy cosmic rays will be affected by the strong 
magnetic fields in the GMZ. Our naive assumptions gloss over 
much of the physics of cosmic ray interaction with molecular 
cloud gas, but they are consistent with typical assumptions 
adopted in the literature and are necessary for a constant cosmic 
ray ionization rate (a Ccr) to be meaningful. 

The effective radius is ofte n computed as the square r oot of 
the total area (in pixels), but |Rosolowsky X Leroy (2006) show 
that this is usually a very bad approximation, so we use the 
Gaussian-corrected (see next footnote) moment-derived radius. 

We correct the linewidth and the radius of dendrogram- 
extracted clumps for clipping assuming the un derlying source 
is Gaussian following Rosolowsky X Blitz (2005 Appendix B). 
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perature To scaled down by an optical depth correction 
1 — exp [—A/'(H 2 )/ 10 ^^ cni“^], and we assume the velocity 
gradient dvjdr = 5 km s~^ pc“^. Of these measurements, 
the linewidth ay is dominant in determining the predicted 
temperature and is the least uncertain. This analysis is per¬ 
formed both for the hand-selected apertures, which have 
arbitrary radii but well-measured linewidths, and the den¬ 
drogram leaf sources. 

The results of the turbulence-predicted temperature 
analysis are shown in Figure which plots our measured 
temperature against the predicted temperature assuming 
turbulent heating is active and the cosmic ray heating is 
negligible {Ccr = 1 x 10“^^ For the hand-extracted 

regions, there is no strong trend in the data, but most data 
points can be explained by turbulent heating alone, per¬ 
haps with cosmic rays setting a high floor temperature. 
There are a few points with measured temperatures below 
the turbulence prediction. Some of these regions include at 
least one velocity component that seems to be genuinely 
cool and broad. A possible explanation for these regions 
is that their linewidths are an indication of bulk motion 
rather than turbulent motion, and therefore they may be 
good candidates for expanding bubble edges, either super¬ 
nova or H ll-region driven. Alternatively, they may contain 
denser gas and therefore be able to cool more efficiently; 
these regions would then be good regions to search for star 
formation. 

The dendrogram-extracted points plotted include only 
those with a Gaussian correction factor < 3, as larger cor¬ 
rections are regarded as unreliable extrapolations. Since 
dendrogram-extracted sources are necessarily coherent in 
position-velocity space, they are less likely than the hand- 
selected features (which were selected in position-position 
space) to represent unassociated line-of-sight features. 

The predicted temperatures for dendrogram clumps fre¬ 
quently exceed the measured temperatures by a large mar¬ 
gin, often a factor of 2 or more. There are two classes of 
explanation: First, that the cooling rate has been underes¬ 
timated. This is entirely plausible, as we have assumed so¬ 
lar neighborhood chemistry (specifically, a CO abundance 
Xco = 10 “^ H 2 , and negligible abundances for other 
coolants), which does not necessarily hold under CMZ con¬ 
ditions. The other possibility is that we have overestimated 
the line width or underestimated the size scale of the clump 
via our extraction method. While we have tried to account 
for this source of error by correcting for clipping, we are 
nonetheless extrapolating our measurements. 

The predicted temperatures for the hand-extracted re¬ 
gions are more uncertain. Most of the hand-extracted re¬ 
gions agree with the prediction at a Icr and all but a few 
at a 2a level. Given that the sizes used for these regions 
are uncertain by a large factor (for square apertures espe¬ 
cially), there is no significant disagreement between Tc^turb 

and Tq^^YY iQdsured' 

To explore the possibility that our cooling rates are un¬ 
derestimated, we examine a particular clump from the cat¬ 
alog, with n = 10^’^ cm“^, a = 8 . 6 /^ = 14.3 km s“^, 
L = 2r = 2(0.588)/^ = 1.65 pc (beam-deconvolved). 
To = 21 K, and To,rad = 1-2 K, where fg is the Gaussian 
correction factor. This clump has a predicted temperature 
Tg = 156 K, but a measured temperature = 91 K. With 
local chemistry, the dominant coolant is CO by 2 orders of 
magnit ude. In DESPOTIC , we s et chemical equilibrium us¬ 
ing the Nelson & Langer (1999) chemical network code and 
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find that the predicted O abundance is increased by an or¬ 
der of magnitude, which in turn makes oxygen a significant 
coolant with Aq ~ fi.bAcoi where A is the cooling rate in 
erg s“^ cm“^ (other coolants remain negligible). The pre¬ 
dicted equilibrium temperature is then Tq = 114 K. While 
this temperature is still a bit high, this example shows that 
the modified chemistry at these high temperatures should 
be accounted for when predicting the gas temperature. 

We may be able to further constrain the heating mech¬ 
anisms using this same data set if more precise velocity 
structure information is obtained for the clouds. Better ve¬ 
locity structure measurements can be achieved by combin¬ 
ing line fitting analysis with systematic clump extraction 
in order to avoid using the extrapolated correction factors. 
This analysis is left for a future work. 


5 . 2 . Temperature variation with spatial scale or radius 


Photon-driven heating mechanisms are subject to attenu¬ 
ation, so clouds with significant X-ray or UV heating may 
have heated shells and cool interiors. Turbulent and cosmic 
ray heating experience no such effects, but they exhibit a 
different dependence on gas density than the corresponding 
cooling mechanisms. We would also expect cooler interiors if 
the clouds are centrally condensed, since more efficient dust 
cooling can dominate heating processes at very high densi¬ 
ties. If clouds are internally heated by gravitational collapse 
and star formation, the cloud interior could be warmer than 
its exterior. 


We have examined the radial profiles of temperature and 
velocity dispersion for selected clouds. We measured the 
ratio Ri and the line width as a function of aperture size by 
extracting progressively larger circular apertures centered 
on the same point. For each spectrum extracted, we have 
fitted the 6-parameter model described in Section 4.2 to 
derive the ratio and linewidth. 


Different variations with scale are observed. For exam¬ 
ple, around ‘The Brick’, the ratio decreases from ^ 0.55 at 
the position i = 0.237°, b = 0.008, which is coincident with 
the highest measured temperature in The Brick (within a 
single spectrum), to ^ 0.33 when averaged over a radius 
^ 70" (Figure pTr). By contrast. Figure 11 3 shows a highly 
turbulent exam^e where the temperature is lowest at the 
center and higher on larger scales. These observations may 
imply that two different heating mechanisms, one internally 
driven (The Brick) and one external (Gl.23-0.08), are ac¬ 
tive within the GMZ. 

These two objects are not meant to be representative of 
the whole CMZ, and the radial properties of clumps need 
to be studied in more detail to draw firm conclusions about 
photon-driven heating processes and thermal balance on 
different size scales. 


5.3. Difference between gas and dust temperatures 


Similar to the results of Ao et al. (2013) for the inner ^ 75 
pc, we find that the P-H 2 CO and dust temperatures are 
not the same anywhere within the inner R < 200 pc. Fig¬ 
ure n!2] shows the fitted P-H 2 CO temperature vs the fitted 


dust temperature from HiGal. The majority of the H 2 CO 
data points are well above the To = Tq line independent 
of the method used to extract the temperature. The few 
sources that are consistent with To = Tq represent an 
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Fig. 9. Plots of gas kinetic temperature vs linewidth. The lines in both plots show modeled temperatures as a function of line 
width given different assumptions as computed by DESPOTIC. A version of this plot with no data is shown in Figure |F.1[ The 
legend shows the various assumptions, with the fiducial assumptions indicated at the top, associated with the dashed gray line; 
these fiducial values are kept fixed except where indicated. The colors of the curves each correspond to a different density, while 
different line styles indicate changes in other physical constants. T he green solid curve shows a pressure equilibrium cas e oc cr“^) 

plotted 


4.2 


with a Larson relation (L oc cr°‘ ) based on|Shetty et al. (| 2012 |. (a) The fitted temperature as described in Section 
against the fitted line width for the by-eye selected regions. The blue symbols are compact ‘clump’ sources and the red symbols 
are large-area ( 8 ^) square regions. The very broad line regions, FWHM > 40 km s“^, are sometimes affected by significant baseline 
issues, and it is possible that the temperature has been underestimated as a result (underestimated 82,1 — 22,0 lines, and therefore 
underestimated temperatures, are possible because of the lower signal-to-noise of the 82,1 — 22,0 line; the 80,3 — 2 o ,2 line is less 
likely to be affected), (b) The same as (a) but for the dendrogram-extracted sources. The points are color coded by signal-to-noise 
in the ratio Ri, with black S/N < 5, blue 5 < S/N < 25, green 25 < S/N < 50, and red S/N > 50. Small, faint points indicate 
ancestors in the dendrogram, while bright, large points indicate leaves. Lower limits on linewidth and temperature are indicated 
with triangles on the edge of the plot. Note that in both panels, any correl ation is even weaker th an that observed in Figure 6 of 


Hiittemeister et al. ( 1998| ), who make the same plot for NH 3 (4,4) and (5,5). Riquelme et al. (2018) also found no correlation using 
NH 3 over a similar range of velocity dispersions. While there is no clear correlation between velocity dispersion and temperature, 
most data points fall within a range that is consistent with turbulent heating if the appropriate size scale is selected; see Figure [Tol 
for plots that account for the differing source sizes. No single cosmic ray ionization rate is capable of reproducing our observations, 
and the higher CRIRs (Ccr > 10“^^ s“^) are inconsistent with most of the data. 


interesting subset of potentially very cool clouds, though 
they all have low signal-to-noise ratios in the present data 
set. Overall, though, our data clearly indicate that, on ^pc 
scales in the CMZ, dust is a coolant and the gas and dust 
are not in equilibrium. This disagreement means t hat gas 
temperatures are inc onsistent with PDR models (Hollen- 
bach fc Tielens||1999 ) and confirms that UV heating is not 
a significant driver of the gas temperatures. 


5.4. Spatial variations in gas temperature and orbit 
comparison 

There are significant spatial variations in the gas tempera¬ 
ture. The most obvious is in and around the Sgr B2 com¬ 
plex, where very high gas temperatures are measured. High 
gas temperatures, Tq ^ 100 — 120 K, were previously mea¬ 
sured in a rid ge to the Galactic nort heast of Sgr B 2 N 

Figure 4b), which 


using CH 3 CN ( |de Vicente et al.||1997 
we now confirm, though our data show Tq ^ 


et al. 


150 K. |OH| 

(|2014[) also sho wed relatively high NH 3 temperatures 


around G859.5, shows cooler temperatures 40K< Tq < 60 
K. Along the ‘dust ridge’ (see Section [h, the high-latitude 
portion is slightly warmer than the lowHatitude region; this 
difference is most evident when comparing The Brick to the 
region just south of it and at I = 0.5°. Figureshows the 
temperature of dendrogram clumps as a function of longi¬ 
tude. The high temperatures {Tq ^ 60 K) persist to high 
longitudes, i ncluding the enigma tic turbulent yet starless 
G1.6 region ( Menten et al.||2009 ). 

In order to search for explanations for the temperature 
variations we observe, we compared our data to predicted 
cloud orbits . Figure shows the orbital fit from |Krmjssen 
et al. ( 2015[ ) and a position-velocity slice through the fitted 


in this hot ridge, and Etxaluze et al. (2018) found a high- 
temperature component (^ 50 — 84 K) in their high-J CO 
spectral line energy distribution analysis. The Brick, Sgr A, 
and the 50 km s“^ and 20 km s“^ clouds exhibit the high¬ 
est temperatures outside of Sgr B2, with Tq > 100 K, in 
ag reement with measure ments in high-J transitions of NH 3 
by Mills & Morris (2018). The western portion of the CMZ, 


temperature cube following that orbit. The orbit does not 
align with all cloud components, especially along the green 
segment, but is still very useful for tracking the properties 
of clouds in common environments. Assuming we can use 
the orbit PPV predictions to locate the clouds in 8 D, we can 
begin to understand the geometry of the CMZ. It appears 
that the far-side clouds (green segment, 1 Myr < t < 2 
Myr), are slightly cooler than many of those on the near 
side, such as “The Brick” (red segment). 

We assigned a time to each dendr ogram-extracted 
clump within 85 km s“^ and 6 pc of the Krnijssen et al. 
(2015) orbit and plotted their temperatures against that 
time in Figure \TE\ We observed a hint of a linear increase 
in temperature with time at the first crossing. The possi- 
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Hand-selected reaions Dendrogram-selected regions 



Turbulence-driven Temperature (K) Turbulence-driven Temperature (K) 

(a) (b) 


Fig. 10. Plots of the p-H 2 CO-measured temperature against the predicted turbulent-driven temperature as discussed in Section 
|5.1| While no trend is obvious in the data, most of the high signal-to-noise data points are fairly close to the To^turb — To,measured 
line. There are few points with To,measured > To^turb, so turbulence generally provides enough - and sometimes too much - energy 
to heat the gas to the observed temperatures. The regions with predicted temperatures higher than observed require additional 
explanation: either their heating rate has been overestimated or their cooling rate underestimated, (a) The hand-selected regions as 
described in Figure the blue symbols are compact ‘clump’ sources and the red symbols are large-area square regions. Triangles 
indicate regions with lower limits on the measured temperature, i.e. any region with a measured 1 — a lower limit on the temperature 
Tg > 150 K. The gray dashed line shows the Tc^turb — To ^measured relation. Regions that are inconsistent with the relation at the 
1 — a level are outlined, (b) The same as (a) for the dendrogram-extracted regions. Regions with a large Gaussian correction factor 
fg>3 are excluded. Regions with Tc^turb > 180 K are indicated with triangles. The points are color coded by signal-to-noise in 
the ratio Ri, with gray S/N < 5, blue 5 < S/N < 25, green 25 < S/N < 50, and red S/N > 50. If the lower signal-to-noise points 
in blue and gray are ignored, the agreement between the predicted temperature and observed temperature is fairly good. 



Fig. 11. (a) Radial plots centered on the southwest portion of ‘The Brick’. The top panel shows the ratio P-H2CO 5'iy(32,i — 

22,o)/*S'iy(3o,3 — 20,2) as a function of aperture size, starting from a single pixel. The derived temperature using Figurej^is shown on 
the right axis; temperatures > 150 K should be regarded as lower limits. The bottom panel shows the line width a — FWHM 12.33 
as a function of aperture size, (b) The same as (a) but for a region centered on Gl. 23 - 0 . 08 , a local peak in the P-H2GO 3o,3 — 2o,2 
emission. In this case, the ratio P-H2GO Su{ 32 ,i — 22,o)/*S'i^(3o,3 — 20,2) is lower at the center and increases toward larger radii. 


ble linear increase in temperature over th e first ~ 0.5 Myr 
would support the Kruijssen et al. (2015) hypothesis that 
the “ring” represents a single gas stream following an ec¬ 
centric orbit with regular pericenter passages by Sgr A*. In 
this model, the pericenter passage drives additional turbu¬ 
lence and induces collapse, which implies increased turbu¬ 
lent energy dissipation (hence heating) after the pericenter 
passage. However, the possible trend is also reasonably con¬ 


sistent with a constant temperature over the same range, 
so we caution against over interpret at ion. 


5.5. Comparison to NH^ measurements from Ott et al. 

(2^ 


Ott et al. (2014) measured the NH 3 ( 1 , 1 ) and (2,2) lines 
over a subset of the region we mapped, from — 0 . 2 ° to 0 . 8 °. 
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Fig. 12. (a) The fitted temperature of region-averaged spectra as described in Section [4^ plotted against the HiGal fitted dust 
temperature. As in Figurethe blue symbols are compact ‘clump’ sources and the red symbols are large-area square regions. The 
black dashed line shows the Tq — Td relation. Nearly all of the data points fall above this relation, and very few are consistent 
with it at the 3cr level, (b) The same as (a) but for the dendrogram-extracted sources, showing a similar discrepancy between dust 
and gas temperature. The points are color coded by signal-to-noise in the ratio i^i, with gray S/N < 5, blue 5 < S/N < 25, green 
25 < S/N < 50, and red S/N > 50. 
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Fig. 13. Derived temperature as a function of longitude colored by LSR velocity for dendrogram-extracted data. The points have 
a size proportional to the 80,3 — 2o,2 line brightness. The more transparent points reflect ancestors in the dendrogram tree, while 
the opaque points show leaves. Lower limits for some of the Sgr B2 regions are shown as red upward-pointing triangles near the 
top of the plot. 


These two lines are a frequently used thermometer sensitive 
to gas temperatures 5 ^ Tq < 40 K, with sensitivity up 
to ^ 80 K but requiring a larger T^ot to Tj^in correction, 
which means the uncertain ties at high temperatures are 
larger ( Mangum et al.|2013 their Figure 1). The NH 3 lines 
have lower critical densities than the H 2 CO 3-2 transitions 
by a few orders of magnitude, ncrit (NH 3 (1,1)) ^ 1 x 10^ 
cm“^ while ricrit(P-H 2 CO 80,3 — 20 , 2 ) ^ 8 x 10 ^ cm“^ and 
UeffiN Rs (1,1)) ^ l e3 while ne//(p-H 2 CO 80,3 - 20 , 2 ) ^ 
The energy levels of the 80,3 — 2 o ,2 


1 X 10^ (|Shirley|2015|f 
and 82,1 — 22,0 lines a 


are 21 and 68 K respectively, while p- 


Ticrit is defined to be the ratio between the spontaneous decay 
coefficient, the Einstein Aui value, to the sum of the collision 
rates into and out of the upper level, rieff is the density required 
to provide a 1 K km s~^ flux density given a temperature and a 


NH 3 (1-1), (2-2), and (4-4) have Eu = 1,42, and 178 K; 
the (4-4) line is therefore required to accurately measure 
temperatures over the range covered by the P-H 2 CO J=3-2 
lines. 


The temperatures derived from the NH 3 line ratio may 
be lower for a variety of reasons. The lower critical density 
means that NH 3 is measuring more diffuse gas (which is 
apparently colder) that is not detected at all in H 2 CO. Since 
H 2 CO only provides measurable temperatures above ^ 50 
K in most locations, there may be low-column gas that 
was too faint to be measured with H 2 CO. However, since 


Glisten et al. (1981) found temperatures comparable to the 
H 2 CO temperatures using higher transitions of NH 3 , there 


reference column density. For either density measure, the species 
with the higher n will trace denser gas. 
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Fi g. 15. The observed te mperature vs time since the most recent pericenter passage of The Brick and the dust ridge clouds on 
the Kruijssen et ah (2015) orbit for dendrogram-extracted clumps within 35 km s~^ and 12 pc of the orbital path, (a) shows the 
entire orbit while (b) shows an expanded view of the first Myr. The colored bar along the bottom of the figure shows the modeled 
orbit velocity at each position. In both panels, the darker symbols represent ‘leaves’ in the dendrogram structure, i.e. compact 
clumps, while the fainter symbols represent their parent structures. There is a hint of an increase in temperature with time, from 
an initial temperature Tg < 50 K to Tg ^ 100 K over a period r ~ 0.8 Myr, as illustrated by the dashed line in panel (b); note 
that the cool features from 0.7-0.9 Myr are at a significantly different velocity and therefore cannot be part of the same stream. 
The Brick is an outlier (it is hotter) whether or not there is a trend. The symbol sizes are proportional to the P-H 2 CO 82,1 — 22,0 
brightness and are meant to give an indication of how reliable the temperature measurements are. The gray vertical lines through 
the circles indicate the formal 1 — a errors on the temperature. 


are either multiple temperature compoueuts iu the diffuse 
gas or peculiar exc itatiou effects affecti ug o ue or both of 

the th ermometers dGiisteu et al. (1981) aud [Morris et al. _ _ 

(1973) suggest that collisioual excitatiou of the ( 2 , 1 ) level work was by Ao et al. ( 2013 ), who observed the iuuer 75 


5.6. Comparison to other measurements with P-H 2 CO 
The largest s urvey of p-H^CO temperatures prior to this 


may de-populate the (2,2) level of NH 3 , for example). The 
pr eseuce of more di f fuse g as traced by at Tgas r\j 350 
K ( Goto et ar]| 2011 [ 2014) suggests that a cold aud diffuse 
compoueut traced by low-J NH 3 is uulikely. 


Figure shows the Ott et al. (2014) temperature map 
iu the same color scheme that we have used iu Figure 
There are some treuds that are commou betweeu the two 
maps: the Sgr A complex aud Sgr B 2 complex are the 
warmest regious, the diffuse gas surrouudiug these is some- 
what cooler, aud the positive latitude compoueut of the 
Moliuari et al.| (|2011|) “riug” is slightly warmer thau the 


uegative latitude compoueut 

The most siguificaut differeuce occurs iu the southwest 
portiou of the Sgr B2 cloud. Iu this regiou, the NH 3 temper¬ 
ature rauges from the high {Tq > SOK) temperatures of the 
Sgr B 2 N hot ridge iuto a cooler backgrouud Tq ^ 30—40 K. 
By coutrast, the H 2 CO temperatures stay high, iudicatiug 
gas temperatures Tq > 150 K exteudiug to the southwest 
of Sgr B 2 S. This differeuce probably iudicates that there 
is a cool, diffuse compoueut beiug detected iu NH 3 , while 
H 2 CO is detectiug a hotter, deuser compoueut. This regiou 
will provide a good comparisou target for other thermomet¬ 
ric observatious. 

Further comparisou of the NH 3 aud H 2 CO thermome¬ 
ters is uecessary, aud will be possible with future ATCA aud 
GBT surveys of higher ammouia trausitious. These will be 
aided b y observatious of both highe r aud lower H 2 CO trau¬ 
sitious (jMaugum & Wootteu 1993) to verify the preseuce 
of a cool, diffuse compoueut. 


pc. They performed a similar aualysis to ours but over 
a smaller area. They reported similar temperatures, with 
some small differeuces due to a differeut set of collisiou 
rates adopt ed (we used 


used 

form. 


Wieseuf eld & Faure 
Ao et a! 


(2013), they 
survey found uni- 


Green (1991)). The 

y high temperatures within the bright Sgr A cloud 
complex. Our survey finds a greater variety in gas temper¬ 
atures due to the larger area covered. 

There have also been ob servatious of p-H^CO temper¬ 


atures with interferometers. | Johnston et al.| (|2014|) noted a 
very high temperature peak, Tq > 320 K, using these lines 
of P-H 2 CO toward The Brick, GO.253+0.015. The ratio of 
the SMA to APEX measurements of both the 80,3 — 2 o ,2 
and 82,1 — 22,0 line in the southwest portion of The Brick is 
0.25, but the ratio Ri ^ 0.45 is approximately the same 
within the r = 25^^ aperture we used. The high tempera¬ 
ture in 
in whic’ 


Johnston et al. (2014) comes from smaller regions 
1 a higher Ri ^ 0.71 is observed. 

Above 300K, the collisiou rates provided by LAMP A 
(Green 1991 Schoier et al.|2005 Wieseufeld fc Faure|2013 ) 
must be extrapolated, which means the models do not pro- 
duce reliable results. Very high Ri values above about 0.7 
cannot be reproduced by high temperatures with the cur¬ 
rent suite of collision rates if the temperature is bounded 
to Tg < 300K, but instead require that the lines be opti¬ 
cally thick and thermalized and therefore that the gas be 
at very high densities. Given the low brightness tempera¬ 
tures observed (~ 2 K even in the SMA ~ 4" beams), if 
the emission is optically thick, it must be coming from an 
extremely small area (it must have a small filling factor). A 
small filling factor for the cloud-scale APEX observations is 
inconsistent with the observed large extent of the emission 
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Fig. 16. (a) A temperature map derived from the NH 3 (l,l)/(2,2) line ratio from Ott et al. (2014 . The NH 3 thermometer has a more limited temperature range when using 
only these two lines, so the map is cut off at Tq — 80 K. (b) The p-H2C0-derived temperature map over the same region from Figure!^ 
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and the low interferometer to sing le-dish ratio. We discuss 
this topic further in Section |5.7[ We conclude, therefore, 
that a lower limit on the gas temperature in The Brick, 
Tg 120 K, is cons i stent w ith our observations, and that 
the Johnston et al. (2014) high Ri measurement should 
be regarded as either a loWr limit Tq ^ 300 K or a sug¬ 
gestion that there is sub-arcsecond optically thick P-H 2 CO 
3o,3 ~ 2o,2 emission at the brightness peaks. 

5.7. Is the gas warm or just optically thick? 


^ are observed, but with much lower filling frac- 
10“^, in The Brick ( Rathborne et al.|2014 


column densities by ^ lOOx. Kauffmann et al. ( |2013 ) ob- 


will dissociate molecules, so in this scenario t he emission 
must come from intermediate velocity shocks (Neufeld & 
Dalgarno 1989). The primary distinction between this in- 


As noted in Section [5^ there is degeneracy in the modeling 
that allows the observed ratios to be produced by very dense 
gas in addition to hot gas. Indeed, ratios approaching R\ = 
1 occur for relatively low temperatures (Tq ^ 30 K) at 
very high densities (n > 10^ cm“^) and very high column 


densities (A^(p-H 2 CO) > 10^^ per km s ^ cm because 
the lines become optically thick and have LTE excitation 
temperatures. 

Since the mean densities at 30" resolution are always 
significantly lower (n < few x 10^ cm“^), such high-density, 
optically thick gas would have to have a very low filling 
factor. The peak P-H 2 CO 3o,3 — 2o,2 brightness ranges from 
0.2 < Tmh ^ 2 K in the detected regions, implying a filling 
factor upper limit // ^ 1 — 7% if we assume the resolved 
brightness temperature equals the dust temperature, Tb = 
Td. In practice, the filling factor would have to be lower 
still, since any optically thin P-H 2 CO in surrounding lower- 
density gas would favor the 3o,3 — 2 o ,2 line and contribute 
to a lower Ri. 

The strongest argument against an optically thick, cold 
gas explanation for our observations is the extremely high 
column density required. At typical P-H 2 CO abundances 
X ^ 10“^, or even more extreme abundances X ^ 10“^, 
the implied local column density is A^(H 2 ) ^ 10^^^^^^ cm“^ 
in a single velocity bin. Column densities up to A/'(H 2 ) 
10^^ cm“ 
tions // - 


termittent model and the more generic turbulent dissipa¬ 
tion heating model is timescale: in the intermittent model, 
the shock-heated gas is not well-mixed with the rest of the 
gas and therefore represents a distinct phase. 

5.8. The CMZ average 

For comparison to extragalactic observations, we have in¬ 
cluded a spectrum averaged over the whole ^ 300 pc ex¬ 
tent of our survey in Figure Eli The peak amplitude is 
50-60 mK and the measured line ratio P-H 2 CO Sj^{32 ,1 — 
‘^2,o)/Su{3o,3 — 20 , 2 ) = 0 . 25 , corresponding to a tempera¬ 
ture Tg' = 65 K at n = 10 ^ cm“^, Tg = 63 K at n = lO^’^ 
cm“^, or Tg = 48 K at n = 10^ cm“^. Over the same area, 
the measured mean dust temperature from HiGal is 23 K, 
though on such large scales it would be more appropriate 
to use Planck or WMAP data to perform this measure¬ 
ment. While the spectral fit does not capture all compo¬ 
nents of the line, the extracted line ratio is representative 



_ _ their 

Figure 4). Even in The Brick, one of the densest clouds in 
the dust ridge besides Sgr B2, higher column densities are 
not observed. So, while it remains possible that optically 
thick P-H 2 CO is responsible for the observed line ratios, it 
is very unlikely, especially outside the densest clouds. 

One remaining possibility is that the P-H 2 CO lines be¬ 
come optically thick because of dramatically reduced tur¬ 
bulent line spreading and therefore increased line trapping. 
Such an effect would require a transition from supersoni¬ 
cally turb ulent to coherent (s nbsonically turbulent) gas mo¬ 
tion (e.g., Pineda et al.|2QlQ ) and could reduce the required 


gas thermal structures in their inner regions. The other 
galaxies in that sample, especially the starburst galaxies, 
have significantly warmer gas and dust, with temperatures 
closer to those measured for Sgr B2 than for the whole 
CMZ. 

M82’s inner regions prese nt an interesting contrast to 
the CMZ.lMiihle et al.M2007j) re port high p-H^CO- derived 
temperature T ^ 200 K, while Weih et al. (2001) found 
NHs-based temperature T ^ 60 K. Taken at face value, 
these measurements imply that both thermometers show 
increasing temperatures in starbursts, but they each trace 
different gas. 


5.9. Implications for star formation 

The star formation properties of the C MZ appear to be dif- 
ferent than in the rest of the Galaxy (|Yusef-Zadeh et M 


2009 2010 Immer et al. 2012 Longmore et al. 2013ap . It 


generally appears defi cient in star formation trac ers relative 


served a few clumps within The Brick with narrow veloc¬ 
ity dispersions that could represent subsonic regions, and 
such clumps could in principle mimic the observed signal 
since their relative line-of-sight velocities are high. How¬ 
ever, there are not enough of these to account for the ob¬ 
served line brightness, and there were also equally many 
high-velocity-dispersion clumps detected. 

A related argument is that the high temperatures come 
from a very small fraction of the gas heated in high veloc¬ 
ity shocks, i.e. that we are seeing interm ittent high tem¬ 
peratures (e.g., Falgarone & Puget 1995). Given the high 
observed velocity dispersions, high velocity shocks capable 
of heating the gas far above the observed temperatures are 
possible. Very high velocity shocks, with Vs > 50km s“^ 


to its dense gas ma ss. Kruijssen et al. (2014) and Krumholz 
& Kruijssen (2015) suggested that one of tne main driving 


factors in this low star formation rate is an increased thresh- 
old for gravitational coll apse caused by turbulent pressure. 
Rathborne et al. (2014) support this claim with observa- 
tions that the gas column probability distribution function 
in The Brick is consistent with a purely turbulent origin. 

Our observations further confirm that turbulence is an 
energetically dominant process in the CMZ. We provide 
an important constraint on the isothermal sound speed. 


We also examine the HOPS ( Walsh et al.pOll ) NH 3 obser¬ 
vations of the GMZ over the same region as our H 2 GO survey, 
and find that the measured (l-l)/( 2 - 2 ) ratio is 1.48, identical 
to that measured for M83 to within measurement errors, so the 
ammonia temperature on these large scales is also consistent. 
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Whole CMZ 



Fig. 17. The spectrum of the whole survey averaged over all pixels. The fit parameters, along with the nominal errors on the 
parameters, are shown in the legend. A single-component fit was used, though many subtler individual components are evident. 
Assuming a constant gas density n = 10' 


,4.5±0.5 


cm the measured ratio implies Tg,cmz = 63l?6 K 


Cs = 0.5 — 0.9 km which is important for setting the 
Mach number that is an essential ingredient in turbulence 
calculations. While the high Cg implies a lower Mach num¬ 
ber, the observed line widths more than make up for it; 
CMZ gas is still more turbulent than Galactic disk gas. 

The high temperature also means that the thermal 
Jeans mass is about a factor of 5 higher than in local clouds 
(assuming Tg,cmz = 60 K, Tq, local = 20 K, Mj^cMz{n = 
10®cm“^) « 9 Mq, Mj^iocaiiji = 10®cm“^) « 2 Mq); 
naively this implies a larger typical fragmentation scal e and 
perhaps a preference toward forming high mass cores ( |Lar-| 
so n|2005 ), though the higher density in CMZ clouds should 


balance this. However, many modern theories of star for¬ 
mation now invoke turbulence to provide the core fragmen- 


for setting core masses ( 

Krumholz & McKee 2005| Hen- 

nebelle & Chabrier||2011l 

Padoan & Nordlund||2011[ |Eeder- 

rath & Klessen|2012 Hennebelle & Chabrier|2013t|Hopkins 


2013). The high temperatures and high Jeans mass provide 


an interesting opportunity to demonstrate that turbulent 
fragmentation, which occurs at the sonic scale, rather than 
Jeans fragmentation is responsible for setting the peak po¬ 
sition of the core mass function. If turbu lence is responsible 


for setting the core fragmentation scale (Offner et al.|2Q13 


Hopkins|2013 ), it should be possible to find sub-Jeans-mass 


in cold, less turbulent local clouds to distinguish the Jeans 
scale from the sonic scale. 

The Equation of State (EOS) in many simulations of 
star-forming regions is generally assumed to follow a fixed 
form t hat matches observations of gas in the solar neighbor¬ 
hood ( Jappsen et al.|[2005 ; Bonnell et al. 2006 Dale et al. 


2012). They assume that Tq = Tp at n ^ 10' 


cm”*^. Our 

measurements suggest that gas at this density is still un¬ 
coupled from the dust, such that the first ‘inflection point’ 
in the EOS should be at a higher characteristic density, 
ric ^ 10^ cm“^. Euture simulations of CMZ-like environ¬ 
ments should account for this difference. 


fragments, and in the warm CMZ it should be easier than 


6. Conclusion 

We present the largest gas temperature map of the CMZ 
to date. We release the data at http: //thedata.harvard, 
edu/dvn/dv/APEX-CMZ- lirnn and the source code for the en¬ 
tire project at https : //github. com/adamginsburg/APEX_ 
CMZ_H2C0, 

The main old result from this paper is that the gas tem¬ 
perature in the galactic center molecular clouds is relatively 
high (50-120 K) and apparently fairly uniform. The temper¬ 
atures whmh__we_d^iye^ higher than the dust tempera¬ 
ture dCiisten et al.||1981l see their conclusion). We have 
added new information on the large scale of this uniformity 


Article number, page 22 of 41 














































Ginsburg et al: APEX CMZ H2CO 


and on the association of these temperatures with particu¬ 
lar physical mechanisms, but overall we confirm those early 
results. There is warm {Tq ^ 60 K) dense gas pervading 
the CMZ out to a radius ^ 200 pc. 

We have examined the thermal balance in CMZ clouds. 
The high but variable gas temperatures throughout the 
CMZ suggest that turbulent heating is the dominant heat¬ 
ing mechanism in Galactic Center dense gas. Cosmic rays 
are not dominant, but may be important in less turbulent 
sub-regions within the CMZ. We are able to place an upper 
limit (^CR < 10 “^^ s“^ in the dense gas because a higher 
CRIR would result in a higher floor temperature than we 
observe. Dust is a coolant in CMZ dense gas, and no other 
heating mechanisms should be able to affect the centers of 
dense clouds. 

There are many regions with expected temperatures 
based on heating from turbulence substantially higher than 
observed. These regions may have chemistry significantly 
different from that in local clouds, with more line coolants 
(e.g., atomic oxygen) available due to the high gas tem¬ 
perature. While we did not investigate the chemistry of 
CMZ clouds in detail, we suggest that exploring chemistry 
in ^ 60 — 100 K, n > 10^ cm“^ material will be essential for 
understanding observations of molecular lines in the CMZ 
and in other galaxies with similar conditions. 

The H 2 CO temperatures are uniformly higher than NH 3 
(l-l)/( 2 - 2 ) temperatures but similar to temperatures de¬ 
rived from higher NH 3 transitions. The lower temperatures 
observed in the lower critical density tracer imply that there 
is a cool, low-density component of the molecular gas. This 
perplexing result is contrary to expectations that the dense 
gas should be the coldest and deserves further study. 

Our new observations are broadly consistent with the 
orbital model of the CMZ dense gas described by |Krmjs-| 


sen et al. (2015). While there are some individual discrep¬ 
ancies, it at least appears that the near-side clouds (The 
Brick and the lettered clouds on the path to Sgr B2) have 
consistent thermal properties. There is a hint that the dust 
ridge clouds fall along a stream of progressively increasing 
temperature as a function of age along the predicted orbits. 
However, the observational uncertainties are also consistent 
with a constant temperature along the ridge, so these data 
cannot be taken as direct proof of the orbital model. 
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Appendix A: Baseline removal 

The baselines in our APEX -1 (SHFI) data were particularly 
problematic, more than is usual in modern heterodyne ober- 
vations. The 218 GHz window we have observed is partic¬ 
ularly sensitive to resonances within the APEX -1 receiver 
that vary on < 1 minute timescales; our off-position calibra¬ 
tions were performed about once per minute and therefore 
were not rapid enough to mitigate this problem completely. 
The baselines can broadly be described as smoothly vary¬ 
ing ripples on the scale of ^ 1/10 the spectral window, 
plus more rapidly varying ripples on 20-40 km s“^ scales. 
In principle, this is a straightforward problem of identifying 
the fourier components associated with each of these scales 
and subtracting them. 

In practice, we discovered that it was not possible to re¬ 
move the dominant baseline structure on either scale with¬ 
out significantly affecting the underlying spectral data. We 
simulated a variety of fourier-space suppression approaches 
by adding synthetic signal to baseline spectra extracted 
from the first few principle component analysis (PGA) com¬ 
ponents of the real spectra. The PGA extraction approach 
is able to pull out the dominant baseline components very 
effectively, but it inevitably includes significant signal in 
the top few most correlated components, especially for the 
strong H 2 CO and lines. We therefore abandoned it 

for the final data reduction. 

For reference, we show an example spectrum that we 
believe to consist entirely of baseline ripples in Figure |AT| 


Appendix B: Spectral fits for regions and apertures 

We provide a catalog of line ratios measured from aperture- 
extracted spectra along with figures showing the best fit 
model for each spectrum. We have extracted spectra in re¬ 
gions with significant P-H 2 CO emission and in commonly 
studied individual regions. This catalog is intended to pro¬ 
vide archival value from our data set and enable future 
studies of individual regions. 

In the hand-extracted source table, three types of region 
are included: 


1 . Circular regions-of-interest (Figure B.li) 

2 . Rectangular regions-of-interest (Figure B.l| 3 ) 

3. 8 ' X 8 ' square observation fields (Fignre jBT] :) 
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Fig. A.l. An example showing some of the worst baselines observed. The plotted spectrum is from an observation on April 2, 
2014, showing the average of the 5% worst spectra. The blue curve shows a a = 5 km s~^ (FWHM= 11.75 km s“^) line centered at 
the 0 km s“^ position of P-H 2 CO 80,3 — 2 o, 2 , illustrating that the baseline ‘ripples’ have widths comparable to the observed lines. 
While we selected the worst 5% in this case, nearly all spectra are affected by these sorts of baselines, and the shape and amplitude 
varies dramatically and unpredict ably. The variation, unpredictable though it is, works in our favor as it averages out over multiple 
independent observations. The 218 GHz region shown here is also the worst-affected; the 220 GHz range that includes the ^^GO 
lines generally exhibits smoother and lower-amplitude baseline spectra. 


The fit tables include multicomponent spectral fits with 
and without spline-based baseline removal (see Section 4.2). 
Spatial parameters for the box (GLON, GLAT, boxheight, and 
boxwidth) or circle (GLON, GLAT, and radius) regions are 
included. The average dust column and dust temperature 
are reported (higaldusttem and higalcolumndens). The 
best fit (maximum likelihood) parameters for gas temper¬ 
ature, P-H 2 CO column, and H 2 density are included along 
with their ±1 — cr marginalized limits. The full column de¬ 
scription is in the tables/README.rst file in the source 
code repository, along with the tables in IPAC format. 

The dendrogram property table is also included in the 
repository and described in the README. 
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Fig. B.l. Figures showing the extracted regions. The background is the integrated masked P-H 2 CO 80,3 — 2o,2 image, (a) Box regions created for large-scale averaging and to 
track the morphology of individual clumps (b) Circular regions selected to highlight points of interest (c) 8' x 8' Box regions matching the observational setup in 2014. 
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Appendix C: Source code and data release 


The reduced and raw data are made available via the 
CfA dataverse (doi:10.7910/DVN/27601). The reduced 
data, including cubes, integrated images, and catalogs, are 
also made available at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/cgi-bin/qcat?J/A+A/ (specific URL to be filled 
in by the Journal). 

The source code for this project in its entirety is avail¬ 
able in the attached tar file and on the internet at https : // 
github.com/adamginsburg/APEX_CMZ_H2C0, Because the 
archives are public, we include scripts to download and pro¬ 
cess the raw data so that all steps of the analysis performed 
here can be performed by any individual with access to 
a computer. We provide a script download_raw_f lies .py 
that will retrieve the raw . apex data from the ESO archive 
and from the dataverse; it requires a valid ESO archive ac¬ 
count to use. The total data reduction process takes ^ 30 
hours on a 48-core, 2012-era linux machine, though most 
steps are not parallelized so the process may be faster on 
more recent machines with fewer cores. About 300 GB of 
free space are required to store the raw, intermediate, and 
reduced data products, though the final products are < 30 
GB, with an actual size depending on whether the baseline- 
subtracted versions are kept separate from the original files. 
The download process may take substantially longer than 
the data reduction. 

To fully reproduce the data products in this paper, some 
effort is required to set up the appropriate directory struc¬ 
ture. The file reduction/run_pipeline_cleese.py pro¬ 
vides a useful template for informing the code about the ap¬ 
propriate directory location, and reduction/README.rst 
includes more detailed instructions. The raw .apex data 
files need to be stored in a common location, and output 
directories need to be specified for the merged data files, 
each epoch’s data (one for the Ao et al. (2013) data, one 
for the 2013 data, and one for the 2014/2015 data), diag¬ 
nostic plots, and the subcubes for each molecule. The func¬ 
tion do_everything in reduction/make_apex_cubes.py 
will then run the pipeline end-to-end. 


for example, very high temperatures require low volume and 
high column densities that occupy a small fraction of the 
total integrated likelihood. 

In the script file example_parplot_constraints.py in 
the associated repository, we examine the systematic un¬ 
certainties in the fitting process. Using a known density, 
temperature, and column, we then use the fitting process 
applied to the dendrogram objects and the hand-extracted 
regions to derive our measured temperatures. Over a range 
10^ < n < 10^ cm“^, the maximum likelihood estimator 
recovers the true temperature to within 25% for all tem¬ 
peratures, and to 10% over the range 10^-^^ < n < 10^-^ 
cm“^. At densities lower than n < 10^ cm“^, the recov¬ 
ered temperature is significantly lower than the true tem¬ 
perature, which is expected since we impose a prior that 
n > 10^ cm“^. Guriously, the expectation value is more bi¬ 
ased than the maximum likelihood estimator, with errors 
reaching 30% and with a much stronger dependence on the 
input temperature. 

Appendix E: Smoothed maps 

In many regions, signal was detected over large areas 
though the peak signal-to-noise ratio remained low. In this 
Appendix, we display smoothed versions of many of the 
main figures in the paper to provide a record of features that 
were weakly detected. The smoothed maps may also give a 
sense of how the temperature varies with spatial scale. 

The cubes were smoothed in the spatial dimensions with 
a (Jfwhm = 33.84" Gaussian to achieve a resolution 45" 
and (Tv,FWHM = 3 km s“^ in the spectral direction. The 
smoothed cubes were also downsampled by a factor of 2 
in the spectral direction (the smoothing and downsampling 
are independent). The same signal extraction analysis de¬ 
scribed in Sect ion was used with a threshold Ta > 3cr. 

Eigure |E.3| shows the peak signal-to-noise in both the 
unsmoothed and the smoothed maps to indicate where the 
temperature measurements can be relied upon. These fig¬ 
ures show the variable-opacity mask used in most of the 
other figures in this paper, with black corresponding to fully 
transparent and white fully opaque. 


Appendix D: A detailed example showing 
parameter constraints 

In Section |4.2[ we discussed the fitting approach used to 
determine a best-fit temperature and uncertainties on that 
temperature. Here, we go through a detailed example show¬ 
ing projections of the full 3D temperature-density-column 
parameter space we explored and the constraints each mea¬ 
surement places in each phase space. 

Eigures pI).H|D.2[ and |D.3| show the 3D parameter space 
of the LVG models viewed from its three unique faces. Each 
figure has 4 panels, each of which shows the effect of one 
of our measurements in constraining the available param¬ 
eter space. All panels include the marginalized likelihood 
contours. Eigure |D^ shows slices through the LVG models 
at the maximum likelihood position. The panels show the 
predicted line brightness of the 3o,3 — 2 o ,2 and 32,i — 22 ,o 
lines in each direction through the LVG parameter cube. 

Eigure [D3] shows the marginalized one-dimensional dis¬ 
tributions for each parameter along with the constraints 
provided by each measurement. Note that the posterior is 
not simply the product of the other displayed distributions; 
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Fig. D.l. The parameter constraints for “cloud c” (Figure^ projected (marginalized) onto the temperature/column density plane, 
(top left) Constraints imposed by the measured ratio p-H^CO Su{^ 2 ,i — ‘^ 2 , 0 )/Su{^o ,3 — 20 , 2 ) are shown in the background filled 
colors, with significance ranges progressing from dark-light from l-3cr. The red line contours show the joint constraints imposed 
by including restrictions on the total column density, volume density, filling factor, and abundance, with 1, 2, 3, and 4a contours 
shown. The x marks the expectation value and the o marks the maximum likelihood value. 

(top right) The same colorscheme as before, showing the constraints imposed by assuming the abundance of P-H 2 CO relative to 
H 2 is as labeled. The abundance does not constrain these parameters, but in Figure [P^ the abundance rules out a substantial 
region of parameter space. 

(bottom left) The same colorscheme as before, showing the constraints imposed by usin g th e me asured mean volume density as a 
lower limit. In this figure, the volume density imposes no constraint, but in Figure [P^ and [P .3 1 it is important. 

(bottom right) The same colorscheme as before, showing the constraints imposed by the measured total column density of H 2 . The 
green contours show constraints imposed by the lower limit on the area filling factor. 
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Fig. D.2. The constraints in density-column parameter space. See Figure [DT] for details. 
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Fig. D.3. The constraints in density-temperature parameter space. See Figure [DT] for details. 
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Fig. D.4. The line brightness of P-H 2 CO 80,3 — 2o,2 (top row) and P-H 2 CO 82,1 — 22,0 (middle row) and P-H 2 CO Su{32,i — 
22,0)/*S'iy(80,3 — 20,2) (bottom row) in the three different projections of parameter space. The grayscale images correspond to a 
slice through the parameter spaces at the location of the b est-fi t parameter. The colored contours show the allowed marginalized 
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Fig. D.5. The one-dimensional probability distributions for each parameter. The legend describes the lines. The blue shaded area 
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Fig. E.l. Same as Figure but smoothed with a 34" FWHM Gaussian. 
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Appendix E.l: Dendrogram cubes 

The dendrogram-extracted catalog was used to build a tem¬ 
perature PPV data cube for visualization purposes. Each 
voxel was assigned a value corresponding to the tempera¬ 
ture in the smallest dendrogram structure it was included 
in. For example, a voxel corresponding to a leaf has the tem¬ 
perature measured for that leaf, while a voxel that is not in 
a leaf but is part of a structure would be assigned the tem¬ 
perature of the smallest structure in which it is included. 
The rest of the cube was assigned NAN values. This cube is 
similar to a smoothed cube, but smoothed over connected 
structures rather than with a symmetric beam. The cube is 
interesting for visualization purposes, but should be treated 
with some skepticism: structures that are connected in PPV 
space in the 80,3 — 2 o ,2 line do not generally have the same, 
constant temperature, but they are forced by construction 
to have the same temperature in this cube. Despite this 



that the dendrogram objects are effectively averages over 
the included area in most cases. Since there is less noise in 
these figures, it is easier to identify interesting features in 
them by eye. 
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Appendix F: Temperature balance plots without 
data 

Figure shows theoretical curves with all data points su¬ 
perposed. For clarity, we re produce that figure with no data 
plotted here in Figure |RH 

The legend shows the various modifications to model pa¬ 
rameters used. The fiducial model is indicated by a dashed 
black line with (cr = 1 x 10“^^ s“^, n = 10^ cm“^, L = 5 
pc, velocity gradient dvjdr = 5 km s“^ pc“^, dust temper¬ 
ature Td = 25 K, dust radiation temperature Tp^rad = 10 
K, and interstellar radiation field Gq. The green solid curve 
shows a model that follows both a size-linewidth relation 
L = pc and is isobaric with n = cm“^, 

where <75 is the velocity dispersion in units of 5 km s“^. 
The green dashed curve shows the same with only the size- 
linewidth relation included. The green dotted curve shows 
a high cosmic-ray ionization rate with no turbulent heat¬ 
ing and is meant to show the absolute floor temperature 
guaranteed by such a high CRIR. 

Appendix G: Source Table 

We include an excerpt from the source tables here. 
The full tables are available in digital form from 
https://raw.githubusercontent.com/adamginsburg/ 
APEX_CMZ_H2C0/vl.l/tables/PPV_H2C0_Temperature_ 
orbit.ipac and https://raw.githubusercontent.com/ 
adamginsburg/APEX_CMZ_H2C0/vl.1/tables/fitted_ 
line_parameters_Chi2Constraints_orbit.ipac, Both 
tables are also available on the dataverse page. A 
complete description of the columns is available at 
https://github.com/adamginsburg/APEX_CMZ_H2C0/ 
blob/vl.1/tables/README.rSt 
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Fig. F.l. Plots of the DESPOTIC-derived equilibrium temperature as a function of line width with no data overlaid. The curves 
are identical to Figure 

Table G.l. H 2 CO Line Parameters and Fit Properties 


Source Name 

Ri 


(Jy 

"^Isr 

Tb{H2CO) 

\og{n{H2)) 

T 

-Lgas 

^gas,turb 





kms“^ 

kms“^ 

K 


K 

K 

Mentenl 

0.35 


0.03 

4.93 ±0.15 

55.85 ±0.15 

1.65 ±0.05 

22.8 

00 rr + ifc).fc)5 

36.8 

Menten? 

0.76 


0.10 

10.34 ±0.97 

54.41 ±0.98 

0.35 ±0.04 

22.4 

SSO.OOtfsf 

92.6 

dn:G0.41+0.05box 

0.29 


0.02 

11.19 ±0.27 

19.59 ±0.27 

1.55 ±0.04 

23.1 

66.90li“« 

102 

20 kms 

0.32 


0.00 

8.15 ±0.02 

17.18 ±0.02 

12.39 ±0.03 

23.3 

fin ^70+12.49 

OZ. ro_20.82 

34.7 

CoolSpot 

0.22 


0.02 

7.26 ±0.39 

27.60 ±0.56 

3.79 ±0.13 

22.9 


140 

Gl.34-0.13 

0.32 


0.02 

5.19 ±0.16 

32.27 ±0.15 

1.68 ±0.05 

22.7 

71X 00+12.49 

41.3 

GO.43-0.05 

0.29 


0.01 

10.25 ±0.13 

89.93 ±0.14 

2.84 ±0.03 

22.5 


87 

G359.49-0.14 

0.25 


0.02 

9.82 ±0.20 

-56.04 ±0.20 

1.57 ±0.03 

22.8 

55.80^9:71 

58.1 

GO.76-0.08 

0.76 


0.05 

12.49 ±0.35 

19.83 ±0.26 

1.62 ±0.10 

23.1 

350.00+f5f 

131 

G0.24-0.05box 

0.20 


0.02 

7.50 ±0.18 

79.72 ±0.18 

1.31 ±0.03 

22.7 

46.08t6;'®l 

56.7 

G0.85-0.04box 

0.31 


0.01 

17.99 ±0.26 

4.03 ±0.27 

1.16 ±0.02 

22.9 

73.84t«3339 

129 

G359.8_box 

0.23 


0.01 

10.53 ±0.25 

11.55 ±0.17 

0.72 ±0.02 

22.5 

53.02t®;®® 

47.5 

G0.83-0.19box 

0.36 


0.01 

10.11 ±0.14 

43.97 ±0.14 

1.76 ±0.03 

22.9 

89.10t^3^j2 

50.3 

Map_002 Off 2 

0.18 


0.04 

11.23 ± 1.33 

6.48 ± 1.54 

0.22 ±0.02 

22.8 


66.6 

Map_007 Off 2 

0.22 


0.01 

26.86 ±0.22 

45.78 ±0.27 

1.30 ±0.01 

23.2 

5l.63t«;33 

142 

Map_012 Off 2 

0.31 


0.06 

9.83 ±1.17 

7.63 ± 1.34 

0.24 ±0.02 

22.5 

79.39tlllt 

59.2 

Map 023 Off 3 

0.08 


0.04 

19.67 ±1.30 

33.04 ± 1.25 

0.28 ±0.01 

22.5 


114 

Map_027 Off 2 

0.15 


0.03 

10.86 ±0.45 

81.72 ±0.43 

0.52 ±0.02 

22.8 

37.76llil 

51.3 

Map_055 Off 3 

0.27 


0.02 

8.36 ±0.27 

10.94 ±0.20 

0.81 ±0.03 

22.8 

62.73+5+ 

47 

Map 116 Off 1 

0.33 


0.03 

13.42 ±0.53 

64.31 ±0.55 

0.41 ±0.02 

22.5 

on 70+16.65 
oU. 1 0 _ 4y 

79.3 

Map 122 Off 1 

0.09 


0.05 

22.05 ±1.61 

10.11 ±1.55 

0.19 ±0.01 

22.7 

30.82^1+® 

128 


Article number, page 39 of 


41 













ASzA proofs: manuscript no. ms 


Table G.2. H 2 CO Parameters and Fit Properties for dendrogram-selected clumps 
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